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ABSTRACT 

SGRJ1550-5418 (previously known as AXP IE 1547.0-5408 or 
PSR J1550-5418) went into three active bursting episodes in 2008 October 
and in 2009 January and March, emitting hundreds of typical Soft Gamma 
Repeater (SGR) bursts in soft gamma-rays. The second episode was especially 
intense, and our untriggered burst search on Fermi/ 'GBM data (8—1000 keV) 
revealed ~450 bursts emitted over 24 hours during the peak of this activity. 
Using the GBM data, we identified a ~150-s-long enhanced persistent emission 
during 2009 January 22 that exhibited intriguing timing and spectral properties: 
(i) clear pulsations up to ~1 10 keV at the spin period of the neutron star 
(P ~ 2.07 s, the fastest of all magnetars), (ii) an additional (to a power- 
law) blackbody component required for the enhanced emission spectra with 
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kT ~ 17 keV, (iii) pulsed fraction that is strongly energy dependent and highest 
in the 50—74 keV energy band. A total isotropic-equivalent energy emitted 
during this enhanced emission is estimated to be 2.9 x 10 40 (D/5 kpc) 2 erg. The 
estimated area of the blackbody emitting region of ~ 0.046(.D/5 kpc) 2 km 2 
(roughly a few xl0~ 5 of the neutron star area) is the smallest "hot spot" 
ever measured for a magnetar and most likely corresponds to the size of 
magnetically-confined plasma near the neutron star surface. 

Subject headings: pulsars: individual (SGR J1550-5418, IE 1547.0-5408, PSR J1550-5418) 
— stars: neutron — X-rays: bursts 



Introduction 



A very small group (roughly half a dozen) of isolated neutron stars have manifested 
themselves in one class as Soft Gamma Repeaters (SGRs) linked by numerous common 
distinguishing properties. Among the most characteristic SGR attributes are (i) X-ray lu- 
minosities much larger (by ~100 times) than the ones expected from their rotational energy 
losses, and (ii) the emission of repeated bursts of soft gamma rays. SGR bursts range from 
"typical" short events lasting ~0.1 s with peak luminosities of L p < 10 41 erg s -1 , to oc- 

and finally 



casional intermediate flares lasting a few seconds with L p ~ 10 42 — 10 43 erg s _1 
to - extremely rare - giant flares lasting a few hundred seconds with L p > 10 45 erg s _1 . 
SGRs were identified together with Anomalous X-ray Pulsars (AXP) as "magnetars": neu- 
tron stars powered by their extremely strong magnet ic fields (surface dipole B ~ 10 14 ~ 15 G; 



(IDuncan fc Thompson 



can 



1992 



Kouveliotou et al. 1998h . Comp r ehens ive reviews on magnetars 



be found in IWoods fc Thompson! fooOfT and iMereghettil (120081 ). and references therein. 



IE 1547.0-5408 was observed with the X-ray Multi-Mirror Mission (XMM-Newton) in 
2004 as a magnetar candidate, selected for its galactic plane location and its relatively 
soft magnetar- like spectrum as seen with the Advance d Satellite for Cosm ology and As- 
trophysics (ASCA) during their Galactic plane survey (jSugizaki et al.l 120011). Although no 
period was detected in the original and follow-up XMM observations, iGelfand fc Gaensler 
(120071 ) also proposed IE 1547.0-5408 as a magnetar candidate based on its spectrum and 
its positional coincidence with an extended galactic radio source G327.24-0.13 (possibly a 
supernova remnant). The subsequent discovery in radio observations of a spin period of 
2.07 s and a period derivative of 2.3 x 10~ n s s _1 led to an estimated dipole surface field of 
B ~ 2.2 x 10 14 G an d confirmed the sou rce's magnetar nature; the source was also renamed 
as PSRJ1550-5418 JCamilo et alil2007h . Its period makes IE 1547.0-5408 the fastest rotat- 
ing magnetar; the source is also one of the only two that emit in radio wavelengths (the other 
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source is an AXP, XTE J1810-197: lHalpern et al.1 12003 ICamilo et alJl2006f ). The distance 



of the source has been estimated by different authors using vari ous methods: ICamilo et a 
( 120071 ) found a distance of ~ 9 kpc by measuring radio dispersion; iGelfand fc Gaenslerl (120071 ) 
estimated ~ 4 kpc assuming an association of the so urce with a possible supernova remnant 



G327.24-0.13; and most recently iTiengo et al.l (120091 ) reported an average distance of 4—5 kpc 
by using observations of an X-ray scattering halo in the Swift/XKT data. Throughout this 
paper we use D 5 = D/5 kpc as the source distance measure. 

On 2008 October 3, IE 1547.0-5408 entered an episode of X-ray activity, emitting several 
typical SGR-like bursts over the next 7 days. During this period, 22 short duration bursts 
were observed with the Gamma-ray Burst Monitor (GBM) on board the Fermi Gamma-ray 
Space Telescope. A detailed analysis of these events is presented in A. von Kienlin et al. 
(2010, in preparation). 

On 2009 January 22, the so urce entered a second period of extremely high X-ray burst 
activity jMereghetti et alJboogh . During the first 24 hours of this "storm" , the Fermi/ GBM 
triggered on the source 41 times: the number of triggers was limited only by the instru- 
ment's capability and did not reflect the actual number of bursts emitted by the source. 
In fact, our on-ground search for untriggered events revealed a total of ~450 bursts during 
this 24 hour period: an unusually high burst frequency from a single source (A.J. van der 
Horst et al. 2010, in preparation). Based on t his SGR-like behaviour, we renamed the source 
as SGRJ1550-5418 kouveliotou et al.ll200gh . 



Upon examination of the data from the first GBM trigger on January 22, we identified 29 
short events riding on an enhancement of the underlying persistent emission lasting ~150 s. 
Closer inspection of this enhancement in different energy ranges revealed periodic oscillations 
with a period consistent with the spin period of SGR J1550-5418. We present here a detailed 
temporal and spectral analysis of this enhanced emission period. In £J5J we briefly describe 
our observations and the GBM instrument and data types. We present our temporal analysis 
results in $31 and our spectral studies in SJU Finally we discuss the physical implications of 
our discovery in §5j 



Instrumentation and Data 



The Fermi/GBM consists of 12 Nal detectors (8—1000 keV) arranged in 4 clusters of 
three each and 2 BGO detectors (0.20—40 MeV) at opposite sides of the spacecraft (for a 
detailed description of the instrument, see iMeegan et al.l 120091 ) . GBM is currently the only 
instrument with continuous broad-band energy coverage (8 keV— 40 MeV) and a wide field 



-4- 



of view (8 sr after taking into account occultation by the Earth) and is, therefore, uniquely 
positioned to accomplish a comprehensive magnetar (or any transient event) monitoring. 
In trigger mode, GBM provides t hree types of data; CTIME Burst, CSPEC Burst, and 



Time Tagged Event (TTE) data flMeegan et all 120091 ) . The CTIME Burst data have a 



time resolution of 64 ms with rather coarse spectral information (8 energy channels). The 
CSPEC Burst data provide high-resolution spectra (128 energy channels) collected every 
1.024 s. Both CTIME Burst and CSPEC Burst accumulate data for ~600 s after a trigger. 
The TTE data provide time-tagged photon event lists for an accumulation time of 330 s, 
starting 30 s prior to the trigger time; this data type provides a superior temporal resolution 
down to 2fis at the same spectral resolution as the CSPEC Burst data. 

The first GBM trigger at the onset of the second active episode from SGR J1550-5418 was 
on 2009 January 22 at 00:53:52.17 UT (= T , GBM trigger number 090122037). In the 600 s 
of the trigger readout we detected many individual short bursts using our on- ground un- 
triggered burst search algorithm. To accept an event as an untriggered burst, we required 
excess count rates of at least 5.5er and 4.5<r in the first and second brightest detectors, respec- 
tively, in the 10—300 keV energy range. We used CTIME data in both continuous (256 ms 
time resolution) and Burst mode (64 ms resolution). Subsequently, we inspected energy- 
resolved burst morphology and compared each detector zenith angle to the source for all 12 
detectors, to determine whether the events originated from SGR J1550-5418. In total we 
identified about a dozen very bright bursts and over 40 less intense bursts within 600 s after 
T (see Figured]). During the same trigger readout we also discovered an enhancement in 
the underlying persistent emission starting at approximately T +70 s and lasting for ~150 s 
(see inset of Figure [I]) . 

One of the events recorded during these 600 s, specifically the burst at T +147 s, was 
so bright that it initiated an Autonomous Repoint Recommendation (ARR), causing the 
spacecraft to start slewing towards the SGR J1550-5418 direction. As the source was already 
close to the boresight of the LAT, the slew angle was pretty small. However, we proceeded to 
check whether the observed emission enhancement was artificially caused by the spacecraft 
slewing. First, we calculated the variation in time of the zenith angle of SGR J1550-5418 for 
each of the 12 GBM detectors. At the onset of the enhancement (T +70 s), the NalO 
detector had the smallest zenith angle to the source of 15°. Due to the ARR, NalO kept a 
constant angle of 18° to the source from T +150 s to T +210 s, after which it constantly 
slewed away from the source until it reached an angle of 23° at To +270 s. During this time 
the persistent emission kept rising until T +150 s, which alone confirms that the enhanced 
emission is intrinsic to SGR J1550-5418. The source was in the field of view of the detector 
until ~To+225 s, at which time it went into an occultation by the LAT. At the same time 
NaI6 was slewing towards the source at an angle of 20°, but the emission was unfortunately 
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obscured by the LAT until ~To+225 s. We note, however, that the enhanced emission was 
not clearly detected with NaI6 after T +225 s. Based on the above, we conclude that the 
rise of the enhanced emission was definitely not caused by the spacecraft slew but by the 
source itself; we cannot unambiguously determine the end of the decay trend (or the total 
duration of the enhancement) in the data due to LAT obscuration. In the analysis presented 
in this paper, we have exclusively used data from Nal (unless noted otherwise), to avoid any 
obscuration effects. Finally, we also checked the LAT data (20 MeV— 300 GeV) of the entire 
day for associated high-energy gamma-ray emission, but found no evidence of high-energy 
photons originating from the direction of SGR J1550-5418. 



3. Temporal Properties of Pulsed Hard X-rays 
3.1. Timing Analysis 



During our search for untriggered events in the first trigger interval of 2009 January 
22 from SGR J1550-5418, we found strong apparent periodic modulations in the enhanced 
emission period from T +130 to 160 s in the 50—102 keV data of detector NalO (see panel 
(c) of Figure [2]). This is the first time to our knowledge that pulsations unrelated to a giant 
flare from a magnetar were clearly seen in the persistent emission of an SGR, in energies 
up to 100 keV. To search for a coherent pulse period, we performed a timing analysis over 
the entire enhancement interval. We first eliminated the times of all the bursts found via 
our untriggered burst search and converted the remaining burst-free intervals to the solar 
system barycenter. For each burst we removed 1 s centered at the burst peak; this elimination 
resulted in a "loss" of ~21 s during the interval T +90— 220 s. As the majority of the SGR 
bursts had durations < 100 ms, our method removed any effect of the burst contributions in 



the t ime series. We then generated a Lomb-Scargle periodogram (]Lomblll975l ; IScargle et al. 



19821 ) over a range of periods from 0.1 s to 10 s using CTIME Burst data in the 50—102 keV 
band. We found a very significant signal with a Lomb power of 72.6 (chance occurrence 
probability, P c ~ 10~ 16 ) at a period of 2.0699 ± 0.0024 s, which is consistent with the spin 
period of SG R J1550-5418. Furthe r, to confirm our detection, we also employed the Z 2 m test 



(with m = 2; iBuccheri et al.lll983l ) on the burst-free and barycentered TTE data. We find a 
coherent signal (with a Z 2 m=2 power of 266, P c ~ 10 -23 ) at the same period. Our spin period 
measurement is consiste nt with the one found for SGR J155 0-5418 using contemporaneous 
X-ray d ata ( Swift /XBTAKuivei et al.ll2009l ; llsrael et al.ll2009l ) and radio data (obtained with 
Parkes; iBurgay et al.l 120091 1 . Therefore, we clearly confirm with the detection of these hard 
X-ray pulsations that the enhanced persistent emission seen in the inset of Figured] originates 
from SGRJ1550-5418. 
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Next we searched in the enhanced persistent emission for evolution in the intensity of 
the pulsations using a sliding boxcar technique. We found that the pulsed signal peaks over 
a 90 s interval, from T +120 to 210 s, which encompasses the peak of the enhancement. 

Finally, we searched for any other intervals exhibiting pulsed emission in the burst-free 
continuous CTIME data of 2009 January 22 and during the four subsequent days, using a 
sliding boxcar of 120 s with 10 s steps. We did not find any additional statistically-significant 
pulsed emission. For the entire search and for all the timing analysis reported here, we used 
more precise spin ephemeris obtained by contemporaneous Swift/XKT, Chandra, XMM- 
Newton and Suzaku observations (G.L. Israel et al. 2010, in preparation). 

3.2. Pulse Profiles 

To investigate the evolution of the pulse profiles with energy, we folded the burst- 
free TTE data spanning 120 s (from T +100 s to T +220 s, which includes the strongest 
pulsation period as found above) with the spin ephemeris of SGR J1550-5418. We estimated 
the background level using the data segment between T to T +60 s. Figure [3] shows the 
source pulse profiles during the enhanced emission interval in six energy bands that have 
the same logarithmic width. The pulse profiles above 110 keV are consistent with random 
fluctuations, and thus not shown. 

Figure [3] indicates that the SGR J1550-5418 pulse profiles in the three lowest energy 
bands are most likely complex (multi-peaked). While the two lowest energy band profiles 
are dominated by the structure around phase 0.7-0.8 (indicated by the dotted lines in Figure 
E]), in the 14—22 keV band we see the emergence of another structure around phase 0.0 
(indicated by the dashed lines in Figure [3]). This pulse becomes equally prominent in the 
22—33 keV range and then dominates in the 33—50 keV band. The pulse profile changes 
remarkably in the 50-74 keV band, which is the most statistically significant of all the energy 
bands investigated, and is distinguished by a broad structure that peaks at around phase 
0.0. The 74—110 keV profile resembles the 50-74 keV one. As noted above, the pulse profile 
above 110 keV is consistent with random fluctuations. Therefore, our results set an observed 
upper energy bound of 110 keV for the hard X-ray pulsations in SGR J1550-5418 during this 
enhanced emission episode. 
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3.3. Pulsed Fraction 

We computed the RMS pulsed fraction using a Fourier based approach as described 



Woods et al.l (120071 ). In summary, we take the Fourier transform of each pulse profile, 
then we calculate the RMS pulsed flux by taking the Fourier coefficients of up to third 
harmonic into account, and finally obtain the pulsed fraction values by dividing the RMS 
pulsed flux by the phase-averaged flux. In Figure IU we show the pulsed fraction spectrum 
of SGR JI550-5418 in the same energy bands as in Figured 

Although marginally significant, there is an indication of a minimum in the RMS pulsed 
fraction around ~30 keV. The RMS reaches its maximum value of 0.55 ± 0.12 in the 50—74 
keV band, and then dips below detection at energies greater than ~110 keV. We will discuss 
the implications of these results along with the results of our spectral analysis in Section £|5j 



3.4. Search for High-Frequency QPOs 



We also searched the period of enhanced emission for any signs of high-frequency quasi- 
periodic oscillations (QPOs) similar to those seen in the aftermath of SGR giant flares and 
attributed t o the excitation of global se ismic modes (llsrael et al.ll2005l ; IStrohmayer &: Watts 
20051 . l200d : IWatts & Strohmaveri l2006h . Using TTE data from the three detectors (NalO, 1 
and 3) with smallest detector zenith angles to the source (and not occulted by the LAT), we 
selected photons with energies less than 100 keV, where the enhanced emission dominates. 
We searched the burst-free data set for periodic and quasi-periodic oscillations (1 Hz and 
2 Hz resolution) and found no significant signals even on timescales as short as 1 s. During 
the period when the emission is strongest (T +100 to 200 s), the 3a upper limit on the 
amplitude of QPOs with frequencies in the range 100—4096 Hz is 7.5% RMS. 

The upper limits are less constraining for frequencies below 100 Hz. A rotational phase 
dependent search also revealed no significant signals. Finally, with the bursts included, we 
searched for short-lived QPOs excited by each event: again, we found nothing significant. 
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4. Spectral Properties of Pulsed Hard X-rays 

4.1. Time-Integrated and Time- Resolved Spectral Analysis 

We analyzed time-integrated and time-resolved spectra of the enhanced emission, using 
the RMFIT (S.lrcl) spectral analysis software developed for the GBM data analysis- 
Similar to the timing analysis, we excluded all bursts identified with the untriggered search 
within the enhancement period. Here, we removed up to 3 seconds of data per burst centered 
at the burst peak, to account for spectral contributions from the wings of each burst; this 
elimination resulted in a loss of ~25 s during the interval T +70— 220 s. We note that 
although some weak bursts may still be included in our enhanced emission spectra, the very 
small intensities of these bursts have practically no effect on our spectral analysis results. 
The spectral study of the untriggered bursts within the enhanced emission indicates that 
their spectra are different from the enhanced emission spectra (A.J. van der Horst et al. 
2010, in preparation). For this analysis, we used only the CSPEC Burst data (8.6—897 keV) 
of detector NalO, which initially had the smallest detector zenith angle to the source (15°) 
and to which the source was visible through most of the enhanced emission. 

Since the Detector Response Matrices (DRMs) of GBM are time dependent due to the 
continuous slewing of the spacecraft, a DRM should be generated for every 2-3 degrees of 
slewing (corresponding to every ~20— 50 s of data). For this analysis, we gen erated DRMs 



for every 50 s starting from To, using GBMRSP vl.7 (see lMeegan et al.l 120091 . for a detailed 
description of the GBM response generation). We used a DRM generated at T + 150 s for 
the time-integrated spectrum (72—248 s), and three DRMs generated at T + 100, 150, and 
200 s, respectively, for the time-resolved spectra: each DRM was centered at the mid-time 
of the accumulation time span of each spectrum. The background spectrum was determined 
by fitting a third-order polynomial function to each energy channel using the burst-free 
intervals (To— 286 to To— 43 s, 1008—1196 s, and 1941—2427 s), with a total accumulation 
time of 896 s. 

We found clear evidence for spectral curvature below 100 keV in the time-integrated 
spectrum of the entire burst-free enhancement period (72—248 s): a single power law thus 
resulted in a very poor fit. We employed five other spectral models; cut-off power law, power 
law + blackbody, optically-thin thermal bremsstrahlung, and single/double blackbody. We 
found that the time-integrated spectrum is best described by a power law + blackbody (see 
Figure [5]). All other spectral models did not provide better fits mainly because they failed 



iR.S. Mallozzi, R.D. Preece, & M.S. Briggs, "RMFIT, A Lightcurve and Spectral Analysis Tool," ©2008 
Robert D. Preece, University of Alabama in Huntsville, 2008 
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to fit the lower energy excess < 10 keV. The best-fit spectral parameters of a power law 
with an additional blackbody are shown in Table [TJ Adding a blackbody (with kT = 18 ± 



4 keV ) to a power law resulted in the most significant improvement in Cash statistics ((Cash 



19791 ) over a single power law (AC-stat = 13.5 for 2 degrees of freedom, corresponding to an 



improvement of 3.25<r). 

The average energy flux over the entire enhancement is (6.5 ± 2.4) x 10 -8 erg cm -2 s _1 
(in 8—150 keV), of which the blackbody component accounts for 19%. As stated earlier, the 
distance estimate to the source is not well constrained; however, assuming a source distance 
of ~ 5 kpc, we estimate a total isotropic emitted energy of 2.9 x 10 40 -D 5 2 erg for the entire 
persistent emission (8—150 keV) during the enhancement. 

To investigate the evolution of the blackbody component and of the source's spectral 
properties in general, we divided the enhanced emission period into three time intervals of 
~50 s each: 74—117 s, 122—169 s, and 173—223 s after the trigger time. The stopping 
time of the last spectrum was T +223 s, because the source was occulted by the LAT for 
NalO around T +225 s. We employed the same set of photon models as the time-integrated 
analysis described above. The first spectrum was best fit by a single power law with no 
evidence of a blackbody or any curvature. The second and third spectra, on the other hand, 
were best described by power law + blackbody models. In the second spectrum (the peak of 
the enhancement) the additional blackbody component was statistically most significant (see 
Figure [5]), and remained significant in the third spectrum as well. The ratio of the blackbody 
flux to the total flux (8—150 keV) was found to be 34% in both intervals. The indices of the 
underlying power-law component, and the blackbody temperature also remained constant, 
at ~— 2.1 and ~17 keV, respectively (within uncertainties; see also Table 1), while the 
power-law amplitude tracked the photon flux. 

We note that the last two time-resolved spectra were also fitted by a cut-off power law 
model with the fit being statistically as good as the power law + blackbody model. The 
difference between the two models becomes apparent only at energies >200 keV, where the 
count rates of our data drop dramatically. For a further comparison, we simulated spectra 
with the best-fit cut-off power law model for the second spectrum, folded them through the 
GBM Nal DRM, and fitted the simulated spectra with a power law + blackbody (and vice 
versa). We found no indication of significant statistical preference between these two models 
due to the low count rates at higher energies. The best-fit cut-off power law parameters for 
the second and third spectra are also shown in Table [TJ 

Finally, we also analyzed the spectrum integrated over the second and third time inter- 
vals (To+122 to 223 s), in which the blackbody component was found to be very significant. 
The blackbody + power law model parameters for this combined spectrum are listed in Ta- 
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ble [H The results were consistent with the time-resolved analysis of the individual spectra 
described above, and the statistical significance of the blackbody component was similar 
to that of the second spectrum (AC-stat = 42.5, corresponding to a 6.2cr improvement). 
A cut-off power law model also provided an adequate fit for this combined spectrum; the 
best-fit parameters of the cut-off power law are also listed in Table HJ 

4.2. Phase-Resolved Spectral Analysis 

We performed spin-phase-resolved spectral analysis of the pulsed enhanced emission, as 
follows: we co-added the burst-free spectrum of each pulse (in T +122 to 223 s, corresponding 
to second and third time-resolved spectra) using TTE data and extracted a phase-maximum 
spectrum (between phases 0.75—1.25 in Figure [3]) and a phase-minimum spectrum (between 
phases 0.25—0.75). The spin phase for each photon was calculated using barycentered times, 
as was done for the timing analysis. We calculated the background spectrum, from the 
burst-free interval at T to T +60 s. 

The spectra of both the phase minimum and maximum were adequately fitted with 
power law + blackbody models, where we kept the power-law indices and the blackbody 
temperatures linked. The values of the linked parameters found in the fit were consistent 
(within la) with those of the time-integrated spectra (see Table [U 122—233 s). However, 
the blackbody component was more significant in the phase-maximum spectrum than in the 
phase-minimum spectrum. The contributions of the blackbody flux to the total flux were 
(52±18)% and (35±18)%, in the phase- maximum and phase-minimum spectra, respectively. 

5. Summary and Discussion 

We report here the discovery of coherent pulsations in the persistent hard X-ray emission 
from SGR J1550-5418 in the Fermi/GBM data lasting ~150 s. Coherent pulsations with 
a 55% RMS pulse fraction have never been detected in the persistent emission at these 
high energies from a magnetar as yet. These pulsations were detected only at the onset of a 
major bursting episode and were not directly related to a major burst or flare from the source. 
The pulse period is consistent with the spin period of SGR J1550-5418 as measured with 
contemporaneous Swift/XRT observations, thus confirming SGR J1550-5418 as the origin of 
the enhanced emission. We estimate the total isotropic-equivalent emitted energy during the 
persistent emission (i.e., excluding burst contributions) to be 2.9 x 10 40 .D 5 2 erg. The thermal 
component accounts for 19% of the total emitted energy; i.e., 5.6 x 10 39 -D 5 2 erg is emitted 
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as blackbody. 



The fact that this enhanced emission was detected at the onset of a major bursting 
episode without evidence of direct association to any particular burst or flare immediately 
before the emission is ver y intriguing. Intermediate flares with p ulsating tails were ob- 
served from SGR 1900+14 (Ibrahi m et al.ll200ll ; iLenters et al.l 120031 ) and very recently from 



SGR J1550-5418 (IMereghetti et al.ll2009i . ~6 hours after T ). Thermal components were also 



found in the decaying tails of i ntermediate events f rom SGR 1900+14 with much lower black- 
body temperatures of ~2 keV (ILenters et al.ll2003l ). The thermal component of the enhanced 
emission we report here is hotter (17 keV), exhibits a strong dependence of the pulse profile 
with energy with a very high RMS pulsed fraction (up to 55%), and is clearly not associated 
with a decaying event tail. Energetically, however, the fluence of this enhanced emission is 
comparable to that of tail emission of the intermediate flares from SGR 1900+14- 

Our timing analysis showed that the detection of pulsations is most significant in the 
120—210 s interval after trigger. We find that the spectrum requires a blackbody component 
along with a power law between 122—223 s, which is consistent with the time interval of the 
most significant detection of pulsations. Moreover, as determined by the energy dependent 
pulse profiles and RMS pulsed fractions, we find that the high-energy pulsations are most 
significant in the 50—74 keV range. Strikingly, the blackbody component of the enhanced 
persistent emission spectrum peaks at around 51 keV (i.e., the Wien peak of 17 keV, see 
Figure [5]). These two independent pieces of evidence lend strong support for a blackbody 
radiation component to account for the curvature in the spectrum of the enhanced emission. 

In our spin-phase-resolved spectral analysis, we find that the blackbody flux to the total 
emission is (52+18) % and (35+18)% in the phase- maximum and phase- minimum spectra, 
respectively. This also suggests that a major contribution to the observed pulsations is from 
the blackbody component. If we assume a surface hot-spot during this pulsating interval, 
then the best-fit blackbody corresponds to an effective radiating area (as projected on the 
plane of the sky, far from the star) of Soo = TcD 2 F/(aT 4 ) fa 0.046-D 5 2 km 2 , where T fa 
2 x 10 8 °K (kT f^ 17 keV) is the observed (gravitationally redshifted) temperature. We 
have used here the blackbody flux at the peak of the pulsations (i.e., phase maximum; 
F fa 5 x 10~ 8 erg cm -2 s _1 ) where the hot-spot is expected to be relatively close to face-on, 
in order to minimize the effects of projection and gravitational lensing by the neutron star, 
so that Soo would be relatively close to the physical area, S, of the hot spot on the neutron 
star surface. For a circular hot-spot, this corresponds to a radius of ~120-D5m. 

The rotational energy of magnetars is insufficient for powering their observed emission, 
since they all have long rotation periods, and their spin-down luminosity is much lower 
than their observed luminosity. Owing to their slow rotation, only a very small fraction 
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(~ i?NS^/c ~ 10~ 5 — 10~ 4 ) of the magnetic f l ux th reading the neutron star corresponds 
to open field lines (IBeloborodov fc Thompson! 120071 ). Since the internal field of magne- 
tars can be significantly stronger and more tangled (or twisted) than the external (largely 
dipole) field, the transfer of magnet ic helicity from the interior to th e exterior of the neu- 
tron star powers magnetar activity (IBeloborodov fc Thompson! 120071 ). As the internal field 
twists t he stellar crust, the magnetosphere als o becomes twisted, possibly in a complex 
manner ([Thompson. Lyutikov. fc Kulkarnil |2002| ) . 



Beloborodov fc Thompson! (120071 ) have shown that the rate of energy dissipation in 
the twisted magnetosphere is L d = J$ e ~ 10 38 Aip(B/10 15 G)(a/.R N s) 2 (e$ e /10 GeV) erg s -1 , 
where I is the net current through the corona, $ e is the voltage along the twisted mag- 
netic lines, a is the size of a twisted region on the stellar surface and Aip characterizes 
the strength of the twist. Identifying a with the inferred size of the hot-spot (~ I2O-D5H1) 
would imply a/R^s ~ 10 -2 which is inconsistent with the observed luminosity of the spot, 
L d ~ 10 38 -D 5 2 erg s" 1 since e<3> e < 10 GeV is expected (limited by pair creation) and Ai[) < 1 
is required for global stability. 

It may be possible for the magnetic twist to grow to a global i nstability level during a 
highl y active bursting period due to frequent starquakes (i.e., Aip > 1; IBeloborodov &: Thompson 
2007). As the magnetosphere untwists, a large amount of energy must be dissipated (ILyutikov 
20061 ). A small "trapped fireball" - plasma of e ± pairs and photons confined by a closed mag- 
netic field region - could then potentially account for the inferred hot-spot, and in particular 
its roughly constant temperature and size. Confining a "fireball" of energy at least compara- 
ble to that emitted by the observed blackbody component, -C^bb 



5.6 x 10 39 £> 5 2 erg, 



ha 3 B 2 > E-, 



iso.bb or B > 



within a region of radius a ~ I2O.D5 m requires Eb(cl) 

1.4 x 10 14 (a/120D 5 m)- 3 / 2 (£ iso> BB/5.6 x 10 39 £> 5 2 erg) 1 / 2 G w 1.4 x 10 14 L> 5 _1/2 G. This is 
consistent with the surface dipole field ofi?~2.2xl0 14 G inferred from the measured PP 
( jCamilo et al.l 120071 ). Therefore, a sufficiently small closed magnetic loop anchored by the 
crust could provide the required confinement. 

Although the neutron star surface is relatively cold, hot spots may naturally form on the 
stellar surface since the energy dissipated in the corona is thermalized as it passes through 
the denser atmosphere and reaches the stellar crust. The large pulsed fraction implies that 
the emitting region responsible for the pulsations, which we identify with one or two hot- 
spots on the stellar surface, is mostly or totally obscured during certain rotation phases. 
This implies that if we are observing emission from two hot-spots (of similar temperature) 
then they cannot be located too far from each other (for example, they cannot be antipodal), 
as most of the stellar surface is visible to an observer at infinity at any given time because 
of strong gravitational lensing by the neutron star. 
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Moreover, we might be observing two hot-spots, possibly corresponding to the two 
footpoints of a twisted magnetic flux tube (or two magnetically confined regions), where the 
second hot-spot is somewhat cooler and dominates the pulsed emission below ~20 keV. This 
might explain the energy dependence of the pulse profiles (see Fig. [3]), with the appearance 
of a second peak at lower energies (below ~33keV, at a phase of ~0.75), as well as the 
increase in the RMS pulsed fraction at the lowest energies (with a local minimum around 
~28keV; see Fig. Hj) and the peak at ~60keV corresponding to the hotter hot-spot. 

This blackbody emission is expected to be accompanied by non-thermal, hi gh-energy 



radiation produced by collisionless dissipation. iBeloborodov fc Thompson! (120071 ) estimate 
that the luminosities in the high-energy and blackbody components should be comparable. 
This is in good agreement with our observations of SGR J1550-5418, where the high-energy 
(power-law) and blackbody contribution to the total luminosity were found to be 65% and 
35%, respectively for the time-resolved spectra. These contributions were 48% and 52% in 
the phase-resolved pulse-maximum spectrum. 

In conclusion, the area of the blackbody emitting region is the smallest "hot spot" mea- 
sured for a magnetar, which likely arises from magnetically confined hot plasma on the neu- 
tron star surfac e, possibly caused by the gradual dissipative untwisting of the magnetosphere 



( lLyutikovll2006l ). If the total radiated energy was initially confined to the inferred extremely 
small size of the enhanced emission region (as in a mini "trapped fireball" scenario), this 
would indicate a very large magnetic energy density (for B > 1.4 x 10 14 £> 5 1/2 G), similar 
to the "trapped fireball" model for the tails of SGR giant flares. The observed enhanced 
emission that we report here is much less energetic than a giant flare tail, while its energy 
is comparable to the tail energy of intermediate events and at the high end of typical SGR 
bursts. Despite some distinct properties, the enhanced emission of SGR J1550-5418 carries 
various flavors of all three SGR phenomena, and thus it is most likely related to the very 
pronounced bursting activity that immediately followed it. 



We would like to thank the referee, Dr. F. Camilo, for his valuable comments. This pub- 
lication is part of the GBM/Magnetar Key Project (NASA grant NNH07ZDA001-GLAST, 
PI: C. Kouveliotou) . We thank G.L. Israel and A. Tiengo for providing the precise spin 
ephemeris and source distance, respectively, prior to their publication. YK and EG acknowl- 
edge EU FP6 Transfer of Knowledge Project "Astrophysics of Neutron Stars" (MTKD-CT- 
2006-042722). JG gratefully acknowledges a Royal Society Wolfson Research Merit Award. 
ER-R thanks the Packard Foundation for support. AJvdH was supported by an appointment 
to the NASA Postdoctoral Program at the MSFC, administered by Oak Ridge Associated 
Universities through a contract with NASA. 



-14- 



REFERENCES 

Beloborodov, A.M. k Thompson, C. 2006, ApJ, 657, 967 

Buccheri, R., et al. 1983, A & A, 128, 245 

Burgay, M., et al. 2009, The Astronomer's Telegram, #1913 

Camilo, F., et al. 2006, Nature, 442, 892 

Camilo, F., et al. 2007 ApJ, 666, L93 

Cash, W. 1979, ApJ, 228, 939 

Duncan, R.C., & Thompson, C. 1992, ApJ, 392, L9 
Gelfand, J.D. & Gaensler, B.M. 2007, ApJ, 667, 1111 
Halpern, J.P., et al. 2005, ApJ, 632, L29 
Ibrahim, A.I., et al. 2001, ApJ, 558, 237 
Israel, G.L., et al. 2005, ApJ, 628, L53 

Israel, G.L., et al. 2009, The Astronomer's Telegram, #1770 
Kouveliotou, C, et al. 1998, Nature, 393, 235 
Kouveliotou, C, et al. 2009, GCN Circular #8915 

Kuiper, L., den Hartog, P.R., Hermsen, W. 2009, The Astronomer's Telegram, #1921 

Lenters, G.T., et al. 2003, ApJ, 587, 761 

Lomb, N.R. 1975, Ap&SS, 39, 447 

Lyutikov, M. 2006, MNRAS, 367, 1594 

Meegan, C, et al. 2009, ApJ, 702, 791 

Mereghetti, S. 2008, Astron Astrophys Rev, 15, 225 

Mereghetti, S., et al. 2009, ApJ, 696 L74 

Scargle, J.D. 1982, ApJ, 263, 835 

Strohmayer, T.E. & Watts, A.L. 2005, ApJ, 632, Llll 



- 15 - 



Strohmayer, T.E. & Watts, A.L. 2006, ApJ, 653, 593 

Sugizaki, M., et al. 2001, ApJS, 134, 77 

Tiengo, A., et al. 2009, ApJ, submitted [arXiv:0911.3"064l yl] 

Thompson, C, Lyutikov, M., & Kulkarni, S.R. 2002, ApJ, 574, 332 

Watts, A.L. & Strohmayer, T.E. 2006, ApJ, 637, L117 

Woods, P.M., & Thompson, C. 2006, in Compact Stellar X-ray Sources, eds. W.H.G. Lewin 
& M. van der Klis, Cambridge Astrophysics Series, 39, p. 547 

Woods, P.M., et al. 2007, ApJ, 654, 470 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 16 - 




j|ig||gj|jgy^^ 



100 200 300 400 500 600 

Time since Trigger (s) 

Fig. 1.— Lightcurve of SGR J1550-5418 in 12-293 keV (GBM NalO CTIME data channels 
1—4). An enlarged view of the pulsed, enhanced emission is shown in the inset. The dashed 
line indicates the background level. 
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Fig. 2. — Lightcurve of SGR J1550-5418 in various energy ranges; (a) 12—27 keV, (b) 
27-50 keV, (c) 50-102 keV, and (d) 102-293 keV. The pulsations are most prominent 
between 50—102 keV (panel c, starting at ~130 s). The bursts have not been removed here 
from the data. The dashed lines in panels (a) to (c) indicate the times of the pulse maxima 
(as calculated using barycentered time). 
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Phase (cycles) Phase (cycles 

Fig. 3. — Pulse profiles of SGR J1550-5418 in equal logarithmic energy intervals; (a) 
10-14 keV, (b) 14-22 keV, (c) 22-33 keV, (d) 33-50 keV, (e) 50-74 keV, and (f) 
74—110 keV. Two cycles are plotted for clarity. The vertical dotted and dashed lines are 
explained in §3.21 
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Fig. 4. — Evolution of RMS pulsed fraction of SGR J1550-5418 as a function of energy. 
Uncertainties are la. The energy bands are the same as those used in Figure [3j 
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Fig. 5. — The photon spectrum of the time interval T +122 to 169 s. The blackbody and 
power-law components are shown separately with dashed curves. The data are binned for 
display purpose only. A 3er upper limit is shown for the highest energy bin. 



- 21 - 



Table 1: Spectral parameters of the enhanced persistent emission period of SGR J1550-5418. 



Time 
since To 


Power Law (PWRL) 1 


Blackbody (BB) 




Energy Flux' 3 

(10~ 8 crgs/cm 2 -s) 




Cut-off Power Law 2 


s 


A 

(X10~ 4 ph/s-cm 2 -kcV) 


7 


N 

(XlO -5 ph/s-cm 2 -keV) 


kT 

(kcV) 


AC-stat 


PWRL BB 


FBB/Ftotal 


Ep Ca k 
(keV) 


72-248 


10.53 (1.96) 


-2.06 (0.10) 


1.23 (0.96) 


17.7 (3.8) 


13.5 


5.30 (2.37) 1.22 (0.28) 


0.19 (0.08) 




74-117 
122-169 
173-223 


5.20 (1.30) 
15.02 (3.05) 
9.18 (3.28) 


-2.15 (0.17) 
-2.09 (0.11) 
-2.14 (0.19) 


No BB 
4.45 (1.58) 
3.49 (2.12) 


17.4 (1.7) 
16.4 (2.7) 


42.9 
15.3 


2.85 (3.30) 

7.82 (4.47) 4.08 (0.65) 
5.05 (3.75) 2.59 (0.72) 


0.34 (0.14) 
0.34 (0.19) 


-1.30 (0.14) 68 ( 7) 
-1.33 (0.25) 59 (10) 


122-223 


13.27 (2.29) 


-2.08 (0.10) 


3.74 (1.50) 


16.5 (1.8) 


42.5 


6.81 (2.99) 2.84 (0.46) 


0.29 (0.10) 


-1.41 (0.13) 65 ( 7) 



1 Power Law Model: /(E)= A(E/100 keV)^ 

2 Cut-off Power Law Model: /(E)= A exp[-E(2 + a)/E pcak ](E/100keV) a 

3 Flux is calculated in 8-150 keV. 

4 1-cr uncertainties are shown in parentheses. AC-stat shows an improvement in Cash statistics 
by adding a blackbody with 2 degrees of freedom to a power law. Cut-off power law parameters 
are shown for the cases where the model also provided adequate fits. 



